Thermodynamic behavior of a phase transition in a model for sympatric speciation 
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We investigate the macroscopic effects of the ingredients that drive the origin of species through 
sympatric speciation. In our model, sympatric speciation is obtained as we tune up the strength of 
competition between individuals with different phenotypes. As a function of this control parameter, 
we can characterize, through the behavior of a macroscopic order parameter, a phase transition 
from a non-speciation to a speciation state of the system. The behavior of the first derivative of 
the order parameter with respect to the control parameter is consistent with a phase transition 
and exhibits a sharp peak at the transition point. For different resources distribution, the transition 
point is shifted, an effect similar to pressure in PVT system. The inverse of the parameter related to 
sexual selection strength behaves like an external field in the system and, as thus, is also a control 
parameter. The macroscopic effects of the biological parameters used in our model reveal thus 
fingerprints typical of thermodynamic quantities in a phase transition of an equilibrium physical 
system. 



The branching of a single population into two or more 
species without prevention of gene flow through geo- 
graphic segregation is known as sympatric speciation 

S S 3 ■ Herbivorous insects have long been consid- 
ered prime candidates for sympatric speciation because 
of an intimate and frequently highly specialized relation- 
ship with their host plants, which serve as habitat, food 
resource, and, often, mating location The apple mag- 
got fly Rhagoletis pomonella has been considered, since 
1966, as the classical example of sympatric speciation in 
progress 'S] . R. pomonella shifted from feeding on the un- 
abscised fruit of its native host hawthorn [Crataegus spp.) 
to utilizing the introduced, domesticated apple [Malus 
pumila) sometime in the mid-1800s in the Hudson River 
Valley region of the state of New York. Genetic evidence 
suggests that the species is in theprocess of shifting and 
adapting to this new host plant |a, |7| ■ 

Two ingredients are important for sympatric specia- 
tion to happen in a population l^ l^lioll: T he competition 
caused by fluctuations in ecology lllUl2l and assortative 
mating caused by selective mating |l3l H^. Ecological 
and sexual selection models have addressed these two as- 
pects of sympatric speciation separately [oj . The starting 
point of ecological models is the assumption that sym- 
patric speciation results from disruptive selection. That 
is, competition for diverse resources leads to separation 
in a population, if individuals with intermediate pheno- 
types are losers when they compete with those with ex- 
treme ones. Such selection can cause sympatric specia- 
tion because it provides an advantage for reproductive 
isolation between opposite, well-adapted, extreme phe- 
notypes, and reproductive isolation can be achieved due 
to evolution of nonrandom mating |l5j . Sympatric spe- 
ciation can also be driven by selective mating, or sexual 
selection, that is, nonrandom mating leading to differen- 



tial reproductive successes of different phenotypes. For 
example if the choice of a mate depends on two traits: 
male display (e.g. nuptial hue, varying from red to blue 
through purple) and female preference for variants of dis- 
play. Some females may prefer red males and others pre- 
fer blue males, this can tear the population apart and 
create a pair of species consisting of red-prefering fe- 
males and red males and of blue-prefering females and 
blue males. 

To study sympatric speciation by simulations we use 
the individual-based Pcnna model 16*. In previous work 
with this model 18, 19, 20] an abrupt ecological change 
was the drive that provoked disruptive selection, as in 
[ll| . which led to speciation through the development 
of assortative mating. A different strategy was used to 
simulate sympatric speciation of predators in a food web 
|2l|. In this case, three types of intra-specific competi- 
tion were adopted, depending on the phenotypic group 
of the predators, and their strength was kept constant 
during the whole simulation. In particular, a parame- 
ter X was introduced, establishing the fraction of the 
populations of extreme phenotypic predators with which 
the intermediate phenotypic individuals would compete, 
besides competing among themselves. In the present pa- 
per we adopt the same kind of constant intra-specific 
competition and study first when speciation is achieved, 
depending on the value of X, for a uniform resource dis- 
tribution per phenotype and a sexual selection of con- 
stant strength. We show that the competition strength 
X plays the role of a control parameter in a phase tran- 
sition, and that the fraction of sexual selective females 
in the population shows behavior similar to an order pa- 
rameter. Furthermore, we show that the transition point 
and its functional form depend on the carrying capacity 
distributions and sexual selection strengths. 
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In the present model, competition for food and assorta- 
tive mating are related to the same phenotypic trait. This 
trait is represented by a new pair of non age-structured 
bit-strings, which are crossed and recombined in the 
breeding process [U- The phenotypic characteristic is 
measured by counting, in this new pair of bit-strings, the 
number of bit positions where both bits are set to 1, plus 
the number of dominant positions (chosen as 16) with 
at least one of the two bits set. It will therefore be a 
number k between and 32, which we will refer to as the 
individual's phenotype. We fix the mutation probability 
per locus, ^ 1, of this phenotypic trait at 0.01. 

In order to consider intra-specific competition depend- 
ing on the individual's phenotype k, we modified the 
logistic Verhulst factor introducing three intra-specific 
competition terms, each one related to a given pheno- 
typic group: 

{Vi{k,t), < fc < ni; especialist 
Vm{k,t), ni < k < n2] intermediate. 
V2{k,t), n2 < k < 32; especialist. 

(1) 

As in the original Penna model, at every time step, and 
for each individual, a random real number uniformly dis- 
tributed between and 1 is generated; if this number is 
smaller than V{k, t), the individual dies. For the extreme 
phenotype groups the competition is given by: 



Vii2)ik,t) 



Pli2)ik,t) + PmiKt) 



F 



(2) 



where Pi(2){k,t) accounts for the population with phe- 
notype k < m (fc > 712) at time t, respectively, Pm{k,t) 
accounts for the population with phenotype fee [rii, ^2], 
and _F is a constant proportional to the carrying capacity, 
taken as 2 x 10^ in our simulations. Individuals with in- 
termediate phenotypes {Pm) compete among themselves 
and also with a fraction X of each population presenting 
an extreme phenotype. The Verhulst factor for them is: 



V,n{k,t) 



Prnik,t)+X *[Pi{k,t)+P2{k,t)] 

F 



(3) 



where X can be thought of as the strength of com- 
petition between intermediate and extreme phenotypic 
populations. Eq[21 means that individuals with extreme 
phenotypes (Pi, P2) compete with those belonging to 
the same phenotypic group and also with the whole in- 
termediate population, but there is no competition be- 
tween extreme phenotypes of different groups because 
we are assuming they are specialized to some extent 
([0,ni — 13), (n2 = 19,32]) on particular resources. 

In order to consider assortative mating, we introduce 
into each female genome a single locus (bit) that codes 
for this selectiveness, also obeying the general rules of 
the Penna model for genetic heritage and mutation. If 
it is set to 0, the female is not selective in mating (pan- 
mictic mating) . It is selective (assortative mating) if this 



locus is set to 1. The mutation probability for this lo- 
cus, which can be in both directions (0 ^ 1), is 0.001. 
Mutated females that are born selective choose mating 
partners according to the following mating strategy: If a 
female has phenotype fc < 16 (fc > 16), it chooses, among 
N„i males, the one with the smallest (largest) phenotype 
value fc; If a selective female has fc = 16 then it chooses 
randomly to act as one of the above. Notice that with 
this strategy all females reproduce every time step from 
age i? = 10 until death. 
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FIG. 1; In (a), (b) and (c) the phenotype distributions of 
the whole population for different X. In the initial steps of 
the simulation, t — 200, the distribution is, in all cases, a 
gaussian centered at intermediate phenotypes. For (a) and 
(b), the distributions at t = 4 x 10^ are equal to those at 
t — ix 10^ and are stationary. For (c) the distribution is not 
stationary neither at 4 x 10^ nor at 4 x 10* (see text), (d) 
The time behavior of the selective females density. 



At the beginning of the simulations females are non- 
selective and all the 6x10'^ individuals (half males and 
half females) have a random phenotype. For the parame- 
ter X = 0, the population with intermediate phenotypes 
does not compete with the extreme phenotypic ones [Eq. 
|3] and, in fact, suffer less competition than the other 
two [Eq. 12]. In this case, nearly all females remain non- 
selective, see Fig^d) lower line, and the phenotype dis- 
tribution corresponds to a stationary gaussian function 
centered at fc = 16, Fig^a) squares. As opposed to the 
situation for X = 0, when we introduce a strong compe- 
tition for the population with intermediate phenotypes, 
by setting, say, X = 1.0, only the individuals with ex- 
treme phenotypes fc = and fc- = 32 survive, as shown 
in Fig^b) squares. In this case the density of selective 
females goes to ps ~ 1 very fast, that is, the females with 
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extreme phenotypes mate only with males of its same 
phenotypic group, see Fig^d) upper line. This means 
that there are two new sympatric species, reproductively 
isolated. For the competition strength X — 0.5, the phe- 
notype distribution is not stationary: in runs that differ 
by the choice of the seed of the random number gener- 
ator, the final distribution sometimes has one maximum 
a,t k = 16, FigOJc) triangles, and some other times it has 
two maxima a.t k — and k — 32, Fig^c) squares. The 
time behavior of the density of selective females presents 
large fluctuations, see Fig^d) central line. Fig. QJi 
shows an important change in the population organiza- 
tion, from a non-speciation state with ps « 0, to a sym- 
patric speciation state with ps 1, as we increase the 
strength of competition, X, for the intermediate pheno- 
types. To determine the range of values of X for which 
sympatric speciation may be obtained, we will analyze 
the behavior of the mean density of selective females, for 
many different strengths of competition. 



the mean size of the whole population is « 50 x 10^, or 
twice the value of the former phase, and the phenotypes 
in the population cluster around only two distinctively 
separated values. Fig. ^b) squares. In the Fig. Efb), 
the large values attained by cr(ps) just above Xc arise 
from large fluctuations in the number of individuals. 

Ecological conditions have been considered an essential 
ingredient for divergence and speciation in sympatry j23l | . 
To evaluate its importance in a phase transition context, 
we simulated different ecological conditions by modifying 
the carrying capacity of the environment, which has so 
far been considered as a constant F in the Eqs. [5] and 
121 It will now be phenotype-dependent and will drive the 
population to experience a disruptive selection between 
the specialist and intermediate phenotypes. Its general 
functional form is: = 2 x 10^ * e"^'^"^^^ Z'^'', where 
each individual, with phenotype k, will feed on a different 
resource niche: For small values of tr^ , the specialists will 
have fewer resources than individuals with intermediate 
phenotype. 
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FIG. 2: a) Mean values the selective females density the or- 
der parameter of the speciation transition, as function of the 
control parameter X. b) Logarithm scale of the standard de- 
viation versus X. For each value of X we have made 10 sim- 
ulations with the same parameters, but using different initial 
seeds for the random number generator. In each simulation 
we calculate the mean value of the density of selective females 



during the last 
of the ten runs. 



10 time steps, and then average the results 
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The behavior of the mean density < > as a function 
of X is shown in FiglS^a). The population changes rather 
abruptly from a non-speciation to a speciation state when 
we change slightly the strength of competition, close to 
Xc = 0.5 Another fingerprint of the macroscopic ef- 
fect, X on < Ps >, is the peak shown by the logarithm of 
the first derivative of the order parameter at Xc, FiglS^b). 
These behaviors are very similar to what happens to an 
order parameter as a function of the control parameter 
in an equilibrium phase transition of a physical system. 
This transition separates a single-species phase from one 
in which two species coexist in sympatry. In the single- 
species phase, Fig. ^a) X < X^ the population presents 
a high diversity with many different phenotypes in the 
population, see Fig^a) squares, and has a mean size of 
« 25 X 10'^. In the two-species phase, Fig. ^a.) X > Xc, 



FIG. 3: The figures show the effect of disruptive selection on 
the speciation transition. The carrying capacity from I to III 
correspond to Fa^ with cr^ — 10^,5 x 10^,10'', respectively. 
Fctj. = F correspond to IV. 

In Fig. Ota) the macroscopic effect of the carrying ca- 
pacity is the shift suffered by the transition point for 
different values of ak, which has then an effect similar 
to pressure in PVT systems. For small value of and 
for X > Xc, case I in Figs. EJa) and (b), the popula- 
tion prefers a non-speciation state, even in the presence 
of a high competition for intermediate phenotypes. This 
happens because there are not enough resources for two 
groups with extreme phenotypes. It is nevertheless im- 
portant to notice that the population has a large diversity 
in this case. That is, the phenotype distribution looks 
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like Fig. ^c) triangles, but it is a stable distribution, see 
Fig. Hb) /. 
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relaxation time since it depends on the Nm value, see 
caption in Figs. 0^ and b. These results, for assortative 
mating are in qualitative agreement with Il4l|. 

We believe these analogies between biological and 
physical parameters will help in the understanding of the 
sympatric speciation process. 

We would like to thank to T.J. P. Penna, J. Nogales, J. 
F. Stilck, E. Curado and P.M.C. de Oliveira for fruitful 
discussions. We acknowledge financial support from the 
agencies CLAF, CNPq, and CAPES. 
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FIG. 4: In (a) and (b) the effect of the sexual selection 
strength, Nm- It is important to point out that the simu- 
lation time for Nm = 3 was 8 x 10^, 20 times bigger than in 
the other cases. 

Sexual selection in the population is associated to the 
number of mating choices each female performs before re- 
production, the parameter N^. The probability of a se- 
lective female with phenotype k < 16 to mate with a male 
of opposite phenotype is P<i6 ~ (0.5)^". In the previ- 
ous section, N^n = 50 and this probability is almost zero, 
meaning that the selective females are highly discrimi- 
natory against the opposite phenotype. With Nm = 3 
the probability becomes P<i6 — 0.125 and it is then pos- 
sible for a selective female with k < 16 to mate with a 
male of phenotype k > 16. When we reduce the num- 
ber of mating choices per female, we can see that the 
phase transition is destroyed. Figs. 01a) and (b). For 
an equilibrium physical system the phase transition dis- 
appears when there is an applied external field, as, for 
an example, happens to the paramagnetic transition of 
magnetic materials at the Curie point. A small value for 
sexual selection strength, equivalent to the application 
of a magnetic field, produces an increase of the selective 
female density, see Fig. EJa) circles for X < X^- This 
density approaches 0.5 since the difference between selec- 
tive mating and panmictic mating is small. 

In conclusion, we reported an investigation of the 
macroscopic effects of the parameters {X competition, 
Fcrj. disruptive natural selection and Nm sexual selection 
strength) on the origin of species by sympatric speciation, 
in a model where fitness and mate choice are represented 
by the same trait. For this model and in the context of 
the phase transition. Fig. El it was possible to quantify 
the ingredients that promote speciation in sympatry. The 
macroscopic effects of the parameters related to natural 
selection by resources are in qualitative agreement with 
other individual-based models that studied the neces sary 
ecological conditions for sympatric speciation [3, IllL Il2l| . 
The effect of sexual selection strength behaves like an ex- 
ternal field in the system and in an physical system a field 
is a control parameter like pressure or temperature. An- 
other characteristic of selective mating is related to the 
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